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Electron-hole Asymmetry and Quantum Critical Point in Hole-doped BaFe2As2 
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We show, from first-principles calculations, that the hole-doped side of FeAs-based compounds is 
different from its electron-doped counterparts. The electron side is characterized as Fermi surface 
nesting, and SDW-to-NM quantum critical point (QCP) is realized by doping. For the hole-doped 
side, on the other hand, orbital-selective partial orbital ordering develops together with checkboard 
antiferromagnetic (AF) ordering without lattice distortion. A unique SDW-to-AF QCP is achieved, 
and J2—J1/2 criteria (in the approximate Ji& J2 model) is satisfied. The observed superconductivity 
is located in the vicinity of QCP for both sides. 

PACS numbers: 74.70.-b, 74.25.Jb, 74.25.Ha, 71.20.-b 
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The superconductivity found in FeAs-layer-based 
compounds [l| is attractive and challenging not only be- 
cause it is the only non-cuprate system, which shows su- 
perconductivity beyond 40 K [2] , but also because of the 
magnetic nature of Fe, and of the multi-band charac- 
ter of the system. Initial studies concentrated on the 
electron-doped compounds, such as i?e(Oi_a;Faj)FeAs [1| 
or i^eOi-^FeAs [3]. For the parent compound LaOFeAs, 
it was first pointed out by Singh et.al [4'] and G. Xu 
et. al. [5], based on first-principles calculations, that 
magnetic instabilities play crucial roles for the under- 
standing of superconductivity, and LaOFeAs is located 
at the border line of both antiferromagnetic (AF) and 
ferromagnetic (FM) instabilities. Soon later it was real- 
ized that significant Fermi surface nesting exists between 
the hole and electron Fermi surfaces (FS) connected by 
a (7=(7r,7r) vector [y, [fl]. Based on transport and op- 
tical measurements and detailed first-principles calcula- 
tions, J. Dong et.al [7| proposed that a stripe-type spin- 
density-wave (SDW) state should be stabilized at low 
temperature, and the essential physics to be discussed 
here is the competing orders between SDW and super- 
conducting states. The predicted SDW state was con- 
firmed by neutron experiment, and the stripe- type spin 
ordering pattern was observed M. The success of LDA 
or GGA-type calculations [9|, |lO| for the present sys- 
tems are in sharp contrast to the failure when applied 
to cuprates, although the electron-correlation may play 
certain roles [ll| . The metallic nature of the SDW state 
also suggests that the FS nesting is an important issue 
to understand the physics. The observed small ordered 
moment [8] (compared with first-principles calculations) 
remains to be a difficult issue, however it was demon- 
strated by recent careful study [9|, [lO, ll2| that the dis- 
crepancy can be much reduced if the optimized structure 
is used. 

The recent interests however move to the hole- 
doped Bai_ajKa^Fe2As2 or Sri_ajKa;Fe2As2 with Tc up to 
38K [13], because very large single crystal can be synthe- 
sized [14]. For the parent compound BaFe2As2, the SDW 
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FIG. 1: The magnetic ordering patterns of two kinds of anti- 
ferromagnetic states. The checkboard ordering is called AFl, 
and the strip-type SDW state [7] is called as AF2 here. Two 
exchange coupUng constants, J\ and J2, are defined for the 
nearest neighbor and next-nearest neighbor interactions, re- 
spectively. 



instability similar to LaOFeAs was observed [l5], and 
it was demonstrated by first-principles calculations ^] 
that BaFe2As2 is electronically similar to LaOFeAs. 
However, as we will show in this paper, the situation 
is true for the parent compounds, but the physics at the 
hole-doping side is quite different with the undoped or 
electron-doped side. 

We performed systematic first principles calculations 
for both hole-doped Bai_a;Ka;Fe2As2 and electron-doped 
LaOi_a;Fa;FeAs, using the virtual crystal approximation 
for the doping. We used the experimental lattice param- 
eters with the optimized internal coordinates of As, as 
suggested by previous studies [l2[ . Two different antifer- 
romagnetic ordered structures (as defined in Fig.l) are 
considered in the present study. 

Parent compound BaFe2As2: Fig. 2 shows the cal- 
culated electronic properties of BaFe2As2 in the non- 
magnetic (NM) solution. Similar to LaOFeAs, there are 
two circle-like hole- type FS around the F point, and two 
ellipse-like electron- type FS around the X point. The X 
point here corresponds to the M point in the Brillouin 
Zone (BZ) of LaOFeAs, because BaFe2As2 has body- 
centered tetragonal structure. On the other hand, dif- 
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FIG. 2: The calculated electronic properties of BaFe2As2: 
(a) The band structure plotted in the Brillouin Zone (BZ) of 
LaOFeAs structure; (b)-(f) The Fermi surfaces (FS), and (f) 
is the plot for total five FS; (g)-(i) The cuted FS for fixed 
kz=0, 7t/2 and tt planes, (j)-(k) The calculated Lindhard 
response function Xo(^) for the qz=0 plane and for the F-X 
line. The X point of BaFe2As2 BZ corresponds to the M point 
of LaOFeAs BZ. 
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FIG. 3: The calculated magnetic moment and stabiliza- 
tion energies of AFl and AF2 states in Bai-ccKccFe2As2 and 
LaOi-ccFa;FeAs, with respect to the NM solution (energy zero 
line) . The lower panels show the ground state phase diagram 
constructed from total energy calculations. Note the two dif- 
ferent quantum critical points (QCP). 



ferent with LaOFeAs, the band dispersion along the z 
direction is stronger. In particular, as shown in Fig. 2 (g)- 
(i), for the kz=0 {kz=7T) plane, the size of hole-pockets is 
smaller (larger) than the electron-pockets, and only for 
the kz=7r/2 plane, the size of two are almost equal. This 
suggests the enhanced three-dimensionality in BaFe2As2. 
The fact that the electron pockets have different orien- 
tations for different kz plane can help us to understand 
the ARPES results, and indeed the calculated FS shape 
can be well compared with the recent ARPES measure- 
ment [17], again suggesting the quality of first-principles 
calculations. 

In addition to the FS characters, we calculated the 
Lindhard response function Xo{q) as shown in Fig. 2 (j)- 
(k). Also similar to LaOFeAs, we found significant FS 
nesting for the q={7r, tt) vector, and the nesting can be 
suppressed by either electron or hole doping. As the re- 
sults, the stripe-type SDW state (called AF2 as shown 
in Fig.l) similar to LaOFeAs is stabilized for the ground 
state (see Fig. 3 for the stabilization energy). The cal- 
culated moment is about 1.5 /i^/Fe, again larger than 
the measure moment of O.S/i^/Fe [15]. All those facts 
suggest that the BaFe2As2 is very similar to LaOFeAs 
except the enhanced three-dimensionality. 

Effect of doping: The electron and hole doped com- 
pounds are however quite different as will be addressed 



in this part. Fig. 3 shows the calculated total ener- 
gies of AFl and AF2 with respect to the NM solu- 
tion, and the corresponding magnetic moment for both 
Bai_a;Ka^Fe2As5 and LaOi-a^F^^FeAs. The qualitative 
behavior as function of doping is definitely very differ- 
ent. We start the discussion from LaOFeAs side. The 
ground state of LaOFeAs is AF2 as suggested before [TJ] 
and confirmed by neutron [8], in the presence of Fermi 
surface nesting. By the electron-doping, the nesting ef- 
fect is suppressed and the stabilization energy of AF2 is 
reduced. It is also shown in Fig. 3 that the size of the or- 
dered moment is reduced by electron-doping. If the dop- 
ing X is approaching 0.2, the AF2 state and its moment 
will be totally suppressed, and the AF2-to-NM quantum 
critical point (QCP) is realized [18]. 

Although the BaFe2As2 itself is quite similar to 
LaOFeAs, however, with the hole-doping, the system be- 
haviors quite differently. First, the energy gains of both 
AFl and AF2 relative to NM state are substantial and 
it is not suppressed significantly by hole-doping. Sec- 
ond, the calculated moment is almost constant and do 
not reduce with doping. What is even more interest- 
ing is that the AFl and AF2 states become degener- 
ate at about x=0.7, and the AFl state is more stable 
than AF2 beyond this doping. The AFl-to-AF2 QCP 
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FIG. 4: Upper panel: The estimated exchange couplings ( Ji 
and J2) as function of hole- doping x in Bai_ccKccFe2As2. Mid- 
dle and lower panels: The calculated magnetic moment for 
each ?>d orbit als for API and AF2 states respectively. 



is therefore reaHzed by hole-doping in sharp contrast to 
the electron-doped side. The constructed ground state 
phase diagram is shown in Fig. 3. The electron and the 
hole-doped sides are asymmetric, and they are physically 
quite different: (1) the AF2-to-NM QCP is realized for 
the electron-doped side, while AF2-to-AFl QCP is real- 
ized for the hole-doped side; (2) the electron-doped side 
can be understood from FS nesting, however the hole- 
doped side behaviors more or less like localized electron. 

The first question to be answered at this stage is why 
the hole-doped side behaviors so differently with the 
electron-doped side? A simple understanding is to look 
at the density of state (DOS) of parent compounds, which 
is asymmetric around Fermi level Ef. It is characterized 
as strong peak and high DOS below Ef^ however dis- 
persive bands (mostly xy) and low DOS above Ef. By 
electron-doping the electron pockets formed by the wide 
xy band dominates, and the system becomes more and 
more itinerant. On the other hand, by hole-doping, the 
Ef is shifted into the narrow band region with high DOS 
peak. The effective Ueff{=U/W) for those flat bands 
are therefore larger than the electron-doped side. In ad- 
dition, we should also notice that the number of d elec- 
trons, which is 6/Fe for undoped compounds, is reduced 
to be more close to 5/Fe (half-filled) by hole-doping. 

The second question is what is the difference between 
the AFl and the AF2 states? The simple answer is that 
orbital-selective partial orbital ordering develops in the 



AFl state. For this purpose, we plot in Fig. 4 the calcu- 
lated magnetic moment decoupled to each orbit als of Fe. 
For the AF2 solution, the magnetic moment is nearly 
uniformly distributed to each orbitals. (The difference 
between yz and zx orbitals is due to the stripe-ordering, 
which breaks the two- fold rotational symmetry). How- 
ever, the AFl solution is unique: one of the orbitals, 
3^2 _ ^2^ which forms a narrow band below Ef^ con- 
tributes to the magnetic moment dominantly. Clearly the 
multi-orbital physics and the orbital degrees of freedom 
play crucial role here. Finally, at the x=1.0, the partial 
orbital ordering of 3z^ — r^ develops in the AFl state. 
With the hole-doping, it is therefore realized that the 
orbital-selective magnetic transition involves (the mag- 
netization is transfered from other orbitals to the 3z^ — r^ 
orbital). 

The third question is why the AFl is stabilized for 
large hole-doping side (relative to AF2 state)? This can 
be understood from the J1-J2 model proposed previously 
for LaOFeAs ^]. For the strong coupling fully- localized 
limit (localized spin system), the J1-J2 model can be used 
to describe the possible ordered states as shown in Fig.l. 
Here Ji is the nearest-neighbor exchange coupling, and 
J2 is the next-nearest-neighbor exchange coupling. The 
positive sign is defined for AF coupling. From the mean 
field solution of this model, the total energy of AFl state 
is given as £^i=-2Ji+2J2, and the energy of AF2 state 
is given as E2=-2J2 per Fe. Therefore if J2=^i/2, the 
two ordered state will be energetically degenerate, and 
strong quantum fluctuation will be expected near this 
critical point. Such model has been suggested to de- 
scribe the physical properties of LaOFeAs, however, first- 
principles calculations [9|, |lO[ suggest that the electron- 
doped LaOFeAs system is actually far away from the 
degeneracy, namely J2 » Ji/2. On the other hand, 
here we show that the quantum degeneracy can be actu- 
ally realized in the hole-doped side. Of course, our sys- 
tem is still away from the fully localized insulating region 
(therefore the usage of J1-J2 model here is not rigorously 
justified). However, We can borrow the simple idea of 
J1-J2 model and qualitative understand the physics. We 
map the calculated total energy to this model, and esti- 
mate the Ji and J2 coupling constants as shown in Fig. 4. 
(The spin S=l is used as suggested from the calculated 
magnetic moment). Indeed, the J2=Ji/2 criteria is real- 
ized at x=0.7, and the stabilization of AFl state as func- 
tion of hole-doping can be understood from the doping- 
dependent modification of Ji and J2 coupling strength. 

Discussions: We would like to address some important 
issues here: 

(1) The proposed AFl solution for large hole-doping 
side does not break the 4-fold rotational symmetry. It 
is therefore expected that the lattice distortion associ- 
ated with the AF2 solution (the SDW state in LaOFeAs, 
which break the 4-fold rotational symmetry) should not 
occur for the AFl phase region. 



(2) In the mean field theory, the quantum phase tran- 
sition between AFl and AF2 happens at Ji = 2J2. For 
a pure 2D system, however, the spacial quantum fluc- 
tuation will completely destroy the long range order in 
the AF2 phase, because of the absence of effective lock- 
ing terms between two different sub-lattices [19|. On the 
other hand, the long range order in AFl state will survive 
under the quantum fluctuation. Therefore beyond mean 
field, the AFl to AF2 transition predicted here is actually 
replaced by AFl to quantum disorder transition. Turning 
on the inter-layer coupling will stabilize the long range 
order for both AFl and AF2 phases in finite tempera- 
ture. Due to the different behaviors of the quantum fluc- 
tuation in AFl and AF2 phases, the AFl phase will be 
more classical with a larger order parameter, while there 
will be still strong quantum fluctuation in the AF2 phase, 
which greatly reduces the order parameter and also may 
stabilize the superconducting phase around the phase 
boundary. This scenario is also supported by the fact 
that the observed ordered moment in BaFe2As2 (0.8//^) 
is larger than that in LaOFeAs (O.S/i^), and closer to 
the calculated one from DFT, because of the enhanced 
3-dimensionality in BaFe2As2. Nevertheless, it should 
be kept in mind that the usage of J1-J2 model is just 
approximate, and also it is not clear yet how the on- 
site multi-orbital physics (inter-orbital fluctuation) may 
modify the story. 

(3) For both the electron and hole-doped sides, the 
observed superconductivity phase region [13] is located 
around the QCP, although two QCPs are physically dif- 
ferent. 

(4) To justify our theory experimentally, it is impor- 
tant to have systematical measurement of spin suscepti- 
bility as function of hole-doping for high-quality single- 
crystals. 

In summary, by first-principle calculation, we show 
that the hole-doped side of FeAs-based compounds is 
very different with its electron-doped side. For the 
high hole-doping side, we predict that the checkboard- 
type AFl state should be stabilized with partial orbital- 
ordering but without lattice distortion. A unique AF2- 
to-AFl QCP is realized with hole-doping, in sharp con- 
trast to the electron-doped side. 

We acknowledge the valuable discussions with Z. Q. 
Wang, J. P. Hu, and the supports from NSF of China and 



that from the 973 program of China (No.2007CB925000). 



[1] Y. Kamihara, et. al., J. Am. Chem. Soc, 

(doi:10.1021/ja800073m). 
[2] X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen and D. 

F. Fang, cond-mat/0803.3603; G. F. Chen, Z. Li, D. Wu, 

G. Li, W. Z. Hu, J. Dong, P. Zheng, J. L. Luo, and N. L. 
Wang, cond-mat/0803.3790; Z. A. Ren, J. Yang, W. Lu, 
W. Yi, G. C. Che, X. L. Dong, L. L. Sun, Z. X. Zhao, 
cond-mat/0803.4283. 

[3] Z. A. Ren, et.al, cond-mat/0804.2582. 

[4] D. J. Singh and M. H. Du, cond-mat/0803.0429. 

[5] G. Xu, W. Ming, Y. Yao, X. Dai, S. C. Zhang, and Z. 

Fang, Europhysics Letters, 82, 67002 (2008). 
[6] LI. Mazin, D.J. Singh, M.D. Johannes, and M.H. Du, 

cond-mat/0803.2740.vl. 
[7] J. Dong, H. J. Zhang, G. Xu, Z. Li, G. Li, W. Z. Hu, D. 
Wu, G. F. Chen, X. Dai, J. L. Luo, Z. Fang, and N. L. 
Wang, Europhysics Letter, 83, 27006 (2008). 
[8] C. de la Cruz, et. al., Nature, 453, 899 (2008). 
[9] T. Yildirim, cond-mat/0804.2252 (2008). 
[10] S. Ishibashi, K. Terakural, and H. Hosono, cond- 
mat/0804.2963 (2008); C. Cao, et.al., cond- 
mat/0803.3236; F. Ma, and Z. Y. Lu, cond- 
F. Ma, Z.Y. Lu, T. Xiang, cond- 
Z. P. Yin, et.al, cond-mat/0804.3355; 
M. A. Korotin, et.al., cond-mat/0804.3453; H. J. Zhang, 
et.al., cond-mat/0803.4487. 
[11] K. Haule, J. H. Shim, and G. Kotliar, cond- 
mat/0803.1279 (2008); K. Haule, G. Kotliar, cond- 
mat/0805.0722; L. Craco, et.al, cond-mat/0805.3636. 
[12] LI. Mazin, and M. D. Johannes, cond-mat/0806.1869. 
[13] M. Rotter, et.al, cond-mat/0805.4630; G. F. Chen, 
et.al., cond-mat/0806.1209; K. Sasmal, et.al., cond- 
mat/0806.1301; G. Wu, et.al., cond-mat/0806.1459; 
[14] G. F. Chen, et.al., cond-mat/0806.2648; J. Q. Yan, et.al., 

cond-mat/0806.2711. 
[15] M. Rotter, et.al, cond-mat/0805.4021; C. Krellner, 
et.al., cond-mat/0806.1043; W. Z. Hu, et.al., cond- 
mat/0806.2652; Q. Huang, et.al., cond-mat/0806.2776. 
[16] I. A. Nekrasov, et.al., cond-mat/0806.2630; F. Ma, et.al, 

cond-mat/0806.3526. 
[17] H. Ding, et.al., cond-mat/0807.0419. 
[18] G. Giovannettil, S. Kumar, J. van den Brink, cond- 

mat/0804.0866. 
[19] C. Fang, H. Yao, W.-F. Tsai, J. Hu, and S. A. Kivelson, 
cond-mat/0804.3843; C. Xu, et.al., cond-mat/0804.4293. 



mat/0803.3286 
mat/0804.3370 



